Implementation of quantum phase gate between two atoms via Rydberg antiblockade and adiabatic passage
Tan Xi, Wu Jin-Lei, Deng Can, Mao Wei-Jian, Wang Hai-Tao, Ji Xin
Department of Physics, College of Science, Yanbian University, Yanji 133002, China

 

† Corresponding author. E-mail: jixin@ybu.edu.cn

Project supported by the National Natural Science Foundation of China (Grant No. 11464046).

Abstract

Combining adiabatic passage and Rydberg antiblockade, we propose a scheme to implement a two-qubit phase gate between two Rydberg atoms. Detuning parameters between frequencies of atomic transitions and those of the corresponding driving lasers are carefully chosen to offset the blockade effect of two Rydberg atoms, so that an effective Hamiltonian, representing a single-photon detuning Λ-type three-level system and concluding the quantum state of two Rydberg atoms excited simultaneously, is obtained. The adiabatic-passage technique, based on the effective Hamiltonian, is adopted to implement a two-atom phase gate by using two time-dependent Rabi frequencies. Numerical simulations indicate that a high-fidelity two-qubit π-phase gate is constructed and its operation time does not have to be controlled accurately. Besides, owing to the long coherence time of the Rydberg state, the phase gate is robust against atomic spontaneous emission.

1. Introduction

Compared with a classical computer, a quantum computer possesses greater capability and faster speed to solve some enormous problems and process massive amounts of information; therefore, it has attracted attention from researchers in recent decades.[13] As is well known, all gate operations in quantum computation can be constructed by a series of elementary one-qubit unitary gates and two-qubit logical gates.[4,5] To put quantum computing into practice, therefore, the quantum gate has been studied theoretically and experimentally in the past several decades.[615] A lot of schemes of implementing quantum gates have been proposed by using diverse physical systems, such as the ion-trap system,[16] linear optical system,[17] cavity-QED system,[1820] NMR system,[21] and superconducting system.[22,23]

As a promising candidate for implementing quantum computation, neutral atoms are of particular interest because of the long coherence time of internal atomic states. Their stable hyperfine energy states are especially suited for encoding logic qubits, and are easily controllable and measurable by using a resonant laser pulse.[24] On the other hand, neutral atoms exhibit large dipole moments when they are excited to Rydberg states, which induces powerful and long-range van der Waals or dipole–dipole interaction. This interaction between the excited Rydberg atoms can lead to Rydberg blockade suppressing resonant optical excitation of over-one Rydberg atoms.[25,26] Many protocols of quantum information processing have been proposed via Rydberg–Rydberg interaction (RRI).[2733] Instead, when a certain relation between Rydberg interaction strength and the detuning between frequencies of atomic transitions and those of corresponding driving lasers is satisfied, the atoms can also be collectively excited to Rydberg states, which results in the so-called Rydberg antiblockade. In the last several years, the Rydberg antiblockade mechanism has been widely applied in quantum information processing and quantum computation.[3443] For example, Su et al. implemented quantum gates in the presence of Rydberg antiblockade;[3638] Shao et al. came up with the Rydberg ground-state antiblockade regime;[39,40] and then Ji et al.[41] and Zhao et al.[42] proposed schemes to fuse entanglement and generate entanglement, respectively, via Rydberg ground-state antiblockade.

In this paper, we propose a scheme for the implementation of a quantum phase gate between two atoms by combining adiabatic passage and Rydberg antiblockade. We carefully design the detuning between frequencies of atomic transitions and those of corresponding driving lasers so as to compensate the blockade effect of two Rydberg atoms. An effective Hamiltonian, which denotes a single-photon detuning Λ-type three-level system and involves the quantum state of two atoms excited simultaneously, is obtained. In the context of the Λ-type three-level effective Hamiltonian, the adiabatic-passage technique is employed to implement a two-qubit phase gate by adopting two time-dependent Rabi frequencies. Adiabatic passage is widely used by slowly varying time-dependent parameters to drive the evolution of a quantum system along its certain eigenstate, generally a dark state with zero eigenenergy, to achieve the desired population transfer. The successful implementation of the two-qubit phase gate by means of adiabatic passage means that the access to the high-fidelity phase gate does not require the precise control of operation time, and that the scheme is robust against atomic spontaneous emission with a very low Rydberg-state dissipative rate.

2. Model and effective dynamics

Figure 1 shows the physical model for implementing the two-qubit phase gate. There are two Rydberg atoms trapped in two separated microscopic dipole traps. Atom 1 possesses one high-lying Rydberg state |r⟩ and two hyperfine ground states |0⟩ and |1⟩, while atom 2 has one more ground state |g⟩ than atom 1. For atom 1, transitions |1⟩1 ↔ |r1 and |0⟩1 ↔ |r1 are driven by classical lasers with Rabi frequencies Ω1 and and corresponding red detuning Δ1 and blue detuning Δ2, respectively. For atom 2, transitions |1⟩1 ↔ |r1 and |g1 ↔ |r1 are driven by classical lasers with Rabi frequencies and Ω2 and corresponding blue detuning Δ2 and red detuning Δ1, respectively. The quantum information is coded in quantum states |0⟩ and |1⟩ of two atoms, and the desired transformation for the target π-phase gate is

Fig. 1. (color online) Sketch of atomic level structures and transitions.

With rotating-wave approximation, the Hamiltonian of the whole system can be written as (assuming ħ = 1):

with Urr being RRI strength dependent on the principal quantum numbers of Rydberg states and the distance between the Rydberg atoms. The evolution of the whole system is governed by Schrödinger equation i|Ψ(t)⟩/∂t = H0|Ψ(t)⟩, in which the state of the whole system at an arbitrary time can be written as the time-dependent superposition of the direct product states of two atoms

in which, for convenience, the time-dependence label “(t)” of the probability amplitude Cmn(t) (m = 0,1,r; n = 0,1,g,r) is dropped. With the direct product states {|mn⟩} being basis vectors and moving H0 in Eq. (2) to a rotation frame with respect to , for which , we have

By imposing the relation Urr = Δ2Δ1 and the large-detuning condition Δ1,2Ω1, , an effective Hamiltonian is induced as[44]

with

in which δk (k = 00, 01, 10, 11, 0g, rr, r0, r1, 1g, 1r, rg, 0r) originates from the Stark shift of the quantum state |k⟩. Ωa(b), , and are effective coupling strengths between |rr⟩ and |11⟩ (|0g⟩), |rg⟩ and |1r⟩, and |0r⟩ and |r1⟩, respectively. By ignoring the terms decoupled to {|00⟩,|01⟩,|10⟩,|11⟩}, the effective Hamiltonian becomes

Then by introducing auxiliary levels and appropriate transitions,[42] the terms of ground-state Stark shifts can be eliminated, and thus the Hamiltonian H2 in Eq. (6) becomes

which is a Hamiltonian of the single-photon detuning Λ-type three-level system. Besides, Hamiltonian Heff involves the quantum state |rr⟩ of two-atom simultaneous excitation, which is caused by the Rydberg-antiblockade effect. Governed by Heff, the state of the effective system at an arbitrary time can be written as

3. Implementing two-qubit π-phase gate via adiabatic passage

Hamiltonian Heff in Eq. (7) does not involve |00⟩, |01⟩, and |10⟩, which indicates that |00⟩ → |00⟩, |01⟩ → |01⟩, and |10⟩ → |10⟩ are accessible readily. If |11⟩ → −|11⟩ is accomplished, the target two-qubit π-phase gate will be implemented. In this section, we next show how to achieve |11⟩ → −|11⟩ by using adiabatic passage.

First of all, we change Ω1 and Ω2 to time-dependent Ω1(t) and Ω2(t), respectively, while other parameters keep time-independent. Hence, Ωa, Ωb, and δrr become time-dependent Ωa(t), Ωb(t), and δrr(t), respectively. The instantaneous eigenstates of Hamiltonian Heff in Eq. (7) are given by

with corresponding eigenvalues being λ+(t) = Ω(t)cotϕ(t), λ0(t) = 0, and λ = −Ω(t)tanϕ(t), respectively, for which we have defined

Based on the theory of adiabatic approximation,[45] slowly varying parameters of the quantum system will constrain the transitions between different eigenstates.[46] If the effective system is in an eigenstate |Φn(t)⟩ (n = +,0,−) initially, then at an arbitrary time the state of the effective system will be[47]

where is the dynamic phase and is the geometric phase. For initial state |11⟩, by setting θ(0) = 0 (assuming the evolution starts at t = 0 and ends at t = tf), the effective system is in |Φ0(t)⟩ initially. It is easy to know that the dynamic phase and geometric phase are both zero, and at an arbitrary time, hence, the state of the effective system is always |Φ0(t)⟩. Therefore, θ(0) = 0 and θ(tf) = π are expected to accomplish the transformation |11⟩ → −|11⟩ (i.e., the target π-phase gate) by following |Φ0(t)⟩.

4. Numerical simulations and discussion

In order to meet the boundary conditions θ(0) = 0 and θ(tf) = π, the Rabi frequencies Ω1(t) and Ω2(t) can be chosen as

in which: T = 0.18tc is the pulse width; τ = 0.12tc is the pulse delay; and tc = 3000/Ω0 is the long pulse time executed to meet the adiabatic approximation. Correspondingly, the effective Rabi frequencies Ωa(t) and Ωa(t) in Eqs. (6) and (7) then can be determined, for which we have chosen other parameters as , Δ1 = 50Ω0, and Δ2 = 10Ω0. With these parameters, we plot the time-dependence of Ω1,2(t) and Ωa,b(t), respectively, in Figs. 2(a) and 2(b). In addition, in Fig. 2(c) we plot the time dependence of θ(t) = arccos[Ωb(t)/Ω(t)]. Obviously, figure 2 (especially Fig. 2(c)) shows that the boundary conditions θ(0) = 0 and θ(tf) = π are satisfied very well.

Fig. 2. (color online) (a) Time-dependence of Ω1(t) and Ω2(t). (b) Time-dependence of Ωa(t) and Ωb(t). (c) Time-dependence of θ(t) = arccos[Ωb(t)/Ω(t)]. Parameters used here: tc = 3000/Ω0, T = 0.18tc, τ = 0.12tc, , , and Δ2 = 10Ω0.

The transformation |11⟩ → −|11⟩ is a key to the implementation of the π-phase gate, and in Fig. 3(a), thus with |11⟩ being the initial state, we give a comparison between the population transfer governed by the total Hamiltonian H0 in Eq. (2) and that governed by the effective Hamiltonian Heff in Eq. (7), for which we have defined Pk = |Ck(t)|2 (k = 11,rr,0g) and . More concretely, in Fig. 3(b), we plot the comparison between the variation of C11(t) governed by H0 and that of governed by the Heff. Figure 3(a) assisted by Fig. 3(b) proves that the evolution governed by Hamiltonian Heff gives the perfect population transfer of implementing the ideal transformation |11⟩ → −|11⟩. Apparently, all curves in Fig. 3 originating from H0, up to slight oscillations, coincide well with those originating from Heff, which declares that the approximation from H0 to Heff is effective and predicts that the implementation of the two-qubit π-phase gate with H0 would be successful.

Fig. 3. (color online) (a) Comparison between the population transfer governed by the total Hamiltonian H0 in Eq. (2) and that governed by the effective Hamiltonian Heff in Eq. (6). (b) Comparison between the variation of C11(t) governed by H0 and that of governed by the Heff. Parameters used here are the same as in Fig. 2.

For a phase gate, the evolution of the phase of a state is a crucial investigated aspect. Therefore, we explore the evolution of the phase for the four different initial states with a complex phase eiϕ0. For the scheme proposed for a π-phase gate, we desire that if the initial state is eiϕ0|00⟩, eiϕ0|01⟩, or eiϕ0|10⟩, neither the state nor the phase will change at all. For the initial state eiϕ0|11⟩, however, after adiabatic passage, the state will keep unchanged but the phase will change with π. In Fig. 4, we give the time-dependence of the phase for the initial state eiϕ0|M⟩ (M = 00,01,10,11) whose evolution is dominated by Hamiltonian Eq. (2), for which we have chosen ϕ0 = π/4 for simplicity. We extract the phase of the state |M⟩ by ϕM = −i ln(CM/|CM|).[48] As shown in Fig. 4, ϕ00, ϕ01, and ϕ10 keep roughly unchanged, i.e., π/4, except slight fluctuations, while ϕ11 is reduced by π with an acceptable approximation after the adiabatic-passage operation, which fits well with expectations. In a word, figure 4 proves that a complex phase involved in the initial state has no influence on the performance of the scheme.

Fig. 4. (color online) Time-dependence of phases for the four different initial states with a complex phase eiπ/4.

We consider a random initial state read as

and after a π-phase gate operation on |Ψ0⟩, the outcome state is

In Fig. 5, we plot the fidelity F(t) = |⟨Ψideal|Ψ(t)⟩|2 by picking some pairs of α and β from α,β ∈ 0 ∼ 2π to test the effectiveness of our scheme with Hamiltonian H0 for implementing the two-qubit π-phase gate. Further more, in Table 1 we list the final fidelity F(tf = 3500Ω−1) with different values of α and β, and all final fidelities in Table 1 are over 0.95. All curves in Fig. 5 and data in Table 1 indicate that the scheme with Hamiltonian H0 for implementing the two-qubit π-phase gate is quite effective.

Fig. 5. (color online) Fidelity F(t) = |⟨Ψideal|Ψ(t)⟩|2 with different pairs of α and β. Parameters used here are the same as in Fig. 2.
Table 1.

Samples of the final fidelity F(tf) with different values of α and β.

.

In fact, a sufficient standard of judging the successful implementation of a quantum gate is the average fidelity defined by

where ρ(t) is the density operator of the quantum system. In the following, therefore, we take atomic spontaneous emission into account, and then consider the average fidelity of implementing the two-qubit π-phase gate by numerically solving the master equation

where , , , , and are Lindblad operators describing the dissipative processes, and γ is the atomic total spontaneous emission rate. In Fig. 6, we plot the time-dependent variation of the average fidelity with different atomic total spontaneous emission rates. As shown in Fig. 6, the greater the atomic total spontaneous emission rate is, the more the final average fidelity descends. Nevertheless, the descent of the final average fidelity is not significant, and even when γ = 0.01Ω0 the final average fidelity is near 0.95. Therefore, we deem that the scheme of implementing the two-qubit π-phase gate is robust against atomic spontaneous emission. After all, the spontaneous emission of Rydberg atoms is far below 0.01Ω0 in practice. In the experiment, the model we proposed to implement the two-qubit π-phase gate can be realized by two 87Rb atoms trapped in two tightly-focused dipole traps.[49,50] Based on the related experimental parameters Ω0 = 2π × 6.8 MHz and γ = 2π × 4.8 kHz ∼ 10−3Ω0, the average fidelity of implementing the two-qubit π-phase gate can reach 98.7%.

Fig. 6. (color online) Time-dependent variation of the average fidelity with different atomic total spontaneous emission rates. Parameters used here are the same as in Fig. 2.
5. Conclusion and prospect

We have accomplished the implementation of the two-qubit phase gate between two Rydberg atoms by combining adiabatic passage and Rydberg antiblockade. An effective single-photon detuning Λ-type three-level Hamiltonian is obtained in the Rydberg antiblockade regime, and a desired two-qubit π-phase gate is able to be executed by using the adiabatic-passage technique. The total Hamiltonian can give a relatively perfect result, which is pretty similar to that given by the effective Hamiltonian. Numerical simulations show that the fidelity with a specified initial state can always be over 95%. The average fidelity with actual experimental parameters is up to 98.7%, which proves that the proposed scheme for implementing the two-qubit phase gate is robust against atomic spontaneous emission.

In the past several years, a lot of methods or techniques, usually called “shortcuts to adiabatic passage” (STAP), have been proposed to speed up an adiabatic passage,[5157] with which many schemes have proposed to achieve quantum information processing and quantum computation.[11,5863] In addition, the combinations between RRI and STAP have been accomplished.[33,40,42] For prospects, therefore, the current scheme is supposed to be accelerated by STAP to construct a fast π-phase gate.

Reference
[1] Deutsch D 1985 Proc. Roy. Soc. London Ser. 400 97
[2] Grover L K 1997 Phys. Rev. Lett. 79 325
[3] Grover L K 1998 Phys. Rev. Lett. 80 4329
[4] DiVincenzo D P 1995 Phys. Rev. 51 1015
[5] Barenco A Bennett C H Cleve R DiVincenzo D P Margolus N Shor P Sleator T Smolin J A Weinfurter H 1995 Phys. Rev. 52 3457
[6] Jaksch D Cirac J I Zoller P Rolston S L Cĉté R Lukin M D 2000 Phys. Rev. Lett. 85 2208
[7] Yi X X Su X H You L 2003 Phys. Rev. Lett. 90 097902
[8] Møller D Madsen L B Mølmer K 2008 Phys. Rev. Lett. 100 170504
[9] Rao D D B Mølmer K 2014 Phys. Rev. 89 030301(R)
[10] Liu X Fang G Y Liao Q H Liu S T 2014 Phys. Rev. 90 062330
[11] Chen Y H Xia Y Chen Q Q Song J 2015 Phys. Rev. 91 012325
[12] Borges H S Villas-Bĉas C J 2016 Phys. Rev. 94 052337
[13] Qin W Wang X Miranowicz A Zhong Z Nori F 2017 Phys. Rev. 96 012315
[14] Duan L M Wang B Kimble H J 2005 Phys. Rev. 72 032333
[15] Deng Z J Zhang X L Wei H Gao K L Feng M 2007 Phys. Rev. 76 044305
[16] Cirac J I Zoller P 1995 Phys. Rev. Lett. 74 4091
[17] Kang Y H Xia Y Lu P M 2016 Quantum Inf. Process. 15 4521
[18] Shao X Q Zhu A D Zhang S Chung J S Yeon K H 2007 Phys. Rev. 75 034307
[19] Zhang Y Q Zhang S 2009 Chin. Phys. 18 4683
[20] Zheng S B 2009 Chin. Phys. 18 3453
[21] Gershenfeld N A Chuang I L 1997 Science 275 350
[22] Yang C P Han S 2006 Phys. Rev. 73 032317
[23] Shao X Q Zhao Y F Zhang S Chen L 2009 Chin. Phys. 18 5161
[24] Saffman M Walker T G Mølmer K 2010 Rev. Mod. Phys. 82 2313
[25] Jaksch D Cirac J I Zoller P Rolston S L Cĉté R Lukin M D 2000 Phys. Rev. Lett. 85 2208
[26] Lukin M D Fleischhauer M Cĉté R Duan L M Jaksch D Cirac J I Zoller P 2001 Phys. Rev. Lett. 87 037901
[27] Shao X Q Zheng T Y Oh C H Zhang S 2014 J. Opt. Soc. Am. 31 827
[28] Shao X Q You J B Zheng T Y Oh C H Zhang S 2014 Phys. Rev. 89 052313
[29] Su S L Guo Q Wang H F Zhang S 2015 Phys. Rev. 92 022328
[30] Shao X Q Wu J H Yi X X 2017 Phys. Rev. 95 022317
[31] Shao X Q Wu J H Yi X X Long G L 2017 Phys. Rev. 96 062315
[32] Zhao P Z Cui X D Xu G F Sjöqvist E Tong D M 2017 Phys. Rev. 96 052316
[33] Kang Y H Chen Y H Shi Z C Huang B H Song J Xia Y 2018 Phys. Rev. 97 042336
[34] Yan D Cui C L Zhang M Wu J H 2011 Phys. Rev. 84 043405
[35] Tian X D Liu Y M Cui C L Wu J H 2015 Phys. Rev. 92 063411
[36] Su S L Liang E Zhang S Wen J J Sun L L Jin Z Zhu A D 2016 Phys. Rev. 93 012306
[37] Su S L Gao Y Liang E Zhang S 2017 Phys. Rev. 92 022319
[38] Su S L Tian Y Shen H Z Zang H Liang E Zhang S 2017 Phys. Rev. 96 042335
[39] Shao X Q Wu J H Yi X X 2017 Phys. Rev. 95 062339
[40] Shao X Q Li D X Ji Y Q Wu J H Yi X X 2017 Phys. Rev. 96 012328
[41] Ji Y Q Dai C M Shao X Q Yi X X 2017 Quantum Inf. Process. 16 259
[42] Zhao Y J Liu B Ji Y Q Tang S Q Shao X Q 2017 Sci. Rep. 7 16489
[43] Chen Y H Shi Z C Song J Xia Y Zheng S B 2018 Phys. Rev. 97 032328
[44] James D F V Jerke J 2007 Can. J. Phys. 85 625
[45] Born M Fock V 1928 Zeitschrift für Physik 51 165
[46] Bergmann K Theuer H Shore B W 1998 Rev. Mod. Phys. 70 1003
[47] Berry M V 1984 Proc. R. Soc. 392 45
[48] Chen Y H Wu Q C Huang B H Song J Xia Y 2016 Sci. Rep. 6 38484
[49] Gaëtan A Miroshnychenko Y Wilk T Chotia A Viteau M Comparat D Pillet P Browaeys A Grangier P 2009 Nat. Phys. 5 115
[50] Miroshnychenko Y Gaëtan A Evellin C Grangier P Comparat D Pillet P Wilk T Browaeys A 2010 Phys. Rev. 82 013405
[51] Chen X Lizuain I Ruschhaupt A Guéry-Odelin D Muga J G 2010 Phys. Rev. Lett. 105 123003
[52] del Campo A 2013 Phys. Rev. Lett. 111 100502
[53] Ibáñez S Chen X Torrontegui E Muga J G Ruschhaupt A 2012 Phys. Rev. Lett. 109 100403
[54] Song X K Ai Q Qiu J Deng F G 2016 Phys. Rev. 93 052324
[55] Chen Y H Xia Y Wu Q C Huang B H Song J 2016 Phys. Rev. 93 052109
[56] Baksic A Ribeiro H Clerk A A 2016 Phys. Rev. Lett. 116 230503
[57] Wu Q C Chen Y H Huang B H Song J Xia Y Zheng S B 2016 Opt. Express 24 22847
[58] Chen Y H Xia Y Chen Q Q Song J 2014 Phys. Rev. 89 033856
[59] Chen Y H Xia Y Chen Q Q Song J 2014 Laser Phys. Lett. 11 115201
[60] Chen Y H Xia Y Song J Chen Q Q 2015 Sci. Rep. 5 15616
[61] Chen Z Chen Y H Song J Xia Y Huang B H 2016 Sci. Rep. 6 22202
[62] Yu L Xu J Wu J L Ji X 2017 Chin. Phys. 26 060306
[63] Xu J Yu L Wu J L Ji X 2017 Chin. Phys. 26 090301